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Multiple  Coulomb  Ordered  Strings  of  Ions 
in  a Storage  Ring 

Rainer  W.  Hasse 


GSI  Darmstadt,  D-6429I  Darmstadt,  Germany 


Abstract.  We  explain  that  the  anomalous  frequency  shifts  of  very  close  masses  measured  in  the 
high  precision  mass  measurement  experiments  in  the  ESR  storage  ring  result  from  the  locking  of 
Coulomb  interacting  strings  of  ions.  Here  two  concentric  strings  which  run  horizontally  close  to 
each  other  for  many  revolutions  are  captured  into  a single  string  if  their  thermal  clouds  overlap. 
They  give  up  their  identity  and  lock  into  an  average  frequency. 


Last  conference  [1],  see  also  ref.  [2],  we  reported  on  theoretical  analyses  of  experi- 
ments with  low  density  highly  charged  ions  cooled  to  ultralow  momentum  spread  in  a 
storage  ring.  It  was  shown  that  the  ions  here  run  on  strings  not  being  able  pass  each  other 
any  more.  In  this  contribution  we  extend  this  work  and  conclude  that  in  recent  high  pre- 
cision mass  measurements  under  similar  conditions  two  (or  more)  adjacent  strings  have 
been  created  which  are  so  closely  situated  that  they  loose  their  identity  and  form  a com- 
mon string,  cf.  also  ref.  [3]. 

In  a recent  paper,  Radon  et  a/.  [4]  reported  on  high  precision  mass  measurements  with 
Schottky  mass  spectrometry  of  many  new  isotopes  in  the  region  57  < Z < 84  with  a high 
resolving  power  of  3.3  x 10^.  This  value  is  not  only  limited  by  the  stability  of  the  bending 
and  focussing  magnets  but  also  by  the  effect  of  anomalous  frequency  shifts  discussed 
below.  In  this  experiment  fully  stripped  isotopes  were  produced  by  fragmentation  of 
a 209gj67+  93Q  MeV^,  separated  in  the  fragment  separator,  and  stored  and 

cooled  in  the  storage  ring  ESR.  The  electron  cooler  cooled  down  the  remaining  few 
thousand  isotopes  to  a longitudinal  velocity  spread  as  low  as  5v/v  = 7 x 10“^,  where 
V = Pc  = 0.685c  is  the  beam  velocity.  The  corresponding  momentum  spread  was  5p/p  — 
Y^5v/v  = 1.3  X 10~^,  where  y—  1.37  is  the  relativistic  factor. 

For  fixed  magnetic  rigidity  Bp  around  7 Tm  each  species  of  the  beam  has  its  own 
frequency  ft  = v/C,-,  where  Q is  the  length  of  its  trajectory  close  to  the  circumference  of 
the  ESR,  C=108.4  m.  Neighboring  masses  Am  by  their  slightly  different  trajectories  due 
to  different  rigidities,  thus  have  different  frequencies 
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Here  q is  the  charge  and  Yt  = 2.67  is  the  transition  energy  which  is  connected  to  the  slip 
factor  Tj  = Y~^  ^ Yt~^  = 0 39  and  to  the  momentum  compaction  factor 
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where  D(z)  is  the  dispersion  function  and  r(z)  is  the  local  radius  of  the  central  orbit. 
Note  that  without  the  effect  of  dispersion  from  the  bending  magnets,  ap  = 0,  the  average 
length  of  the  trajectories  and,  hence,  all  frequencies  would  be  the  same. 

The  measurements  were  done  at  the  16th  harmonic  of  the  revolution  frequency 
1.9  MHz,  hence  f=30.4  MHz  and  the  noise  frequency  from  the  last  term  of  eq.  (1)  is 
about  15  Hz.  At  present,  experiments  are  under  way  [5]  at  the  transition  energy,  y = Yi- 
in  which  these  noise  effects  are  eliminated  although  no  cooling  is  needed.  Each  spec- 
trum obtained  was  averaged  up  to  ten  thousand  times  and  the  masses  were  identified  and 
determined  with  sophisticated  analysis  on  the  basis  of  eq.  (1). 

In  the  results  typical  deviations  were  found  from  the  measured  frequencies  of  known 
isotopes  as  compared  to  the  known  calibration  values.  As  shown  in  Fig.  1,  if  there  arc 
at  least  two  masses  so  close  to  each  other  that  their  difference  in  frequencies  is  smaller 
than  about  80.. .90  Hz,  on  the  average  the  lower  (higher)  one  is  shifted  characteristically 
to  a higher  (lower)  value.  According  to  eq.  (1)  this  anomalous  effect  limits  the  mass 
resolution  to  about  20  fm  and  such  masses  then  have  to  be  discarded  from  the  results. 
They  are  the  major  reason  for  the  systematic  errors.  The  authors  believe  that  due  to 
a high  number  of  ions  circulating  with  close  distances  next  to  each  other  in  the  ESR, 
Coulomb  interactions  may  lead  to  collective  coupled  motion. 


frequency  difference  [Hz] 


FIGURE  1.  Average  anomalous  frequency  shifts  of  neighboring  isotopes  vs.  difference  between  their 
measured  frequencies.  Each  point  is  the  average  of  about  100  data  (after  ref.  [4]).  The  result  of  the  model 
calculation  included  (full  line)  and  the  sawtooth  curve  (dashed  line,  only  partly  in  the  frame)  which  would 
result  without  folding  are  discussed  below. 

In  what  follows  we  explain  and  make  a model  for  the  origin  of  these  anomalous  fre- 
quency shifts  on  the  basis  of  the  anticipation  that  the  ions  run  on  correlated  trajectories. 
The  situation  here  is  similar  to  the  one  of  the  machine  experiments  on  the  anomalous 
jump  to  very  low  momentum  spreads  in  the  ESR  [6]  and  in  the  SIS  [7j.  There  it  has  been 
shown  [2]  that  under  the  experimental  conditions  intra-beam  scattering  is  completely 
suppressed  and  that  the  ions  cannot  pass  each  other  any  more  and,  thus,  run  on  strings. 

With  the  parameters  give  above  and  a mean  tune  of  Q = 2.3  the  Wigner-Seitz  radius 
is  flws  = (3ir//2m(Op)’^^  = 14  g/n,  where  the  betatron  frequency  C0|3  = QcOrev  is  the 
Q-fold  of  the  revolution  frequency.  The  Wigner-Seitz  radius  roughly  is  the  closest 
distance  which  ions  can  attain  in  infinite  Coulomb  matter  at  zero  temperature.  With 
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2000  (5000)  ions  in  the  ring  the  particles  have  an  average  distance  in  the  beam  direction 
of  J = 5(2)  cm  and  a linear  density  of  A.  = a^s/d  = 0.00025  (0.0007),  which  is  well 
within  the  region  of  stable  strings,  X < 0.7  [8].  The  average  Coulomb  energy  between 
two  particles  at  distance  d is  ec~  q^y/d=0A9  (0.47)  meV.  With  the  methods  of  ref.  [2] 
we  derive  a longitudinal  kinetic  energy  Tjong  = /n(c|3  5v/ v)^/ (8  In  2)=3  meV  and  take  the 
same  transverse  kinetic  energy  as  in  [6],  rtrans=l-5  eV.  The  FWHM  transverse  thermal 
diameter  of  the  projection  of  the  ion  cloud  onto  the  transverse  plane  calculated  in 
the  harmonic  oscillator  potential  generated  by  the  betatron  frequency  (Op  = 2nQ^c/C 

becomes  25x  = 4^/^/A,^^trans/3cc  = 46  ptm  [9].  An  estimate  for  the  transverse  emittance 
from  etrans  = (§.^)^/Ptrans  with  Ptrans  « 15  m [6]  being  the  average  transverse  beta 
function  in  the  free  sections  yields  the  extremely  low  value  4 x 10“^  n mm  mrad.  All 
these  parameters  are  of  the  same  order  of  magnitude  as  in  the  experiments  of  ref.  [6]. 
Consequently,  in  the  present  experiment  the  ions  have  to  run  on  strings  as  well. 

Nuclides  with  different  masses  but  with  the  same  velocity  run  on  different  trajectories 
which  are  horizontally  displaced  so  that,  apart  fom  fluctuations  from  noise,  v = Cf 
remains  constant,  thus 
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This  difference  in  length  of  trajectory  transforms  into  the  horizontal  displacement 
-Ax  = (C/2n)A///  = 17  m X A///  which  attains  the  value  43. ..48  /tm  for  the  experi- 
mental value  of  Fig.  1 of  the  onset  of  the  anomalous  frequency  shifts  at  A/=80...90  Hz. 
In  addition,  the  dispersion  by  the  bending  magnets  also  causes  a horizontal  displacement 
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FIGURE  2.  Schematic  drawing  of  the  projection  of  the  ion  clouds  onto  the  transverse  plane  in  the 
unlocked  and  locked  cases. 

due  to  the  remaining  thermal  fluctuations.  With  the  dispersion  function  of  the  ESR, 
the  average  dispersion  is  about  D « 2.4  m.  Hence  the  maximum  thermal  displacement 
Dbv/v  ss  1.7/rm  should  be  small  as  compared  to  the  displacement  of  eq.  (3).  An  addi- 
tional displacement  from  the  centrifugal  force  can  be  neglected  since  this  effect  would 
be  smaller  by  a factor  of  Ax/C  » 4 x 10~^. 

This  fact  leads  to  the  model  of  Fig.  2,  namely  that  if  the  displacement  is  larger  than  the 
thermal  diameter,  the  two  (or  more)  strings  run  on  well  separated  trajectories  yielding 
two  (or  more)  distinct  peaks  in  the  Schottky  spectrum  which  lie  at  the  correct  positions. 
However,  in  the  opposite  case  when  the  clouds  overlap,  the  strings  loose  their  identity 
and  lock  partially  or  fully  into  the  same  frequency,  thus  smaller  frequencies  apparently 
become  larger  and  vice  versa. 
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FIGURE  3.  Theoretical  averaged  Fourier  spectra  around  the  third  harmonic  from  Molecular  Dynamics 
calculations  with  30  particles  in  a repeating  boundary  box  distributed  over  two  strings  with  the  occupation 
numbers  written  at  the  bases.  Unlocked,  partially  locked  and  locked  means  that  the  displacements  of  the 
strings  are  larger,  in  between  and  smaller  than  the  thermal  diameter,  ©long  and  0,rans  are  the  longitudinal 
and  transverse  kinetic  energies  in  units  of  the  Coulomb  interaction  energy  of  two  particles  at  distance  d 
and  X is  the  linear  density. 
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Assuming  that  this  locking  is  not  ad  hoc  but  by  virtue  of  the  averaging  procedure  over 
the  experimental  data  that  there  is  a smooth  transition  with  a probability  proportional  to 
the  overlap  region  of  the  clouds,  P = [a  — sin(a)]/7t,  where  a is  the  opening  angle  of 
Fig.  2,  and  folding  this  probability  together  with  the  experimental  resolution  of  15  Hz 
into  the  sawtooth  curve  — /0(|9O  Hz|  — /)  of  shifts  (the  dashed  lines  in  Fig.  1),  where 
0 is  the  Heaviside  function,  yields  the  result  of  Fig.  1.  The  dashed  sawtooth  line  would 
result  if  the  strings  were  always  captured  if  their  thermal  radii  start  to  overlap.  Note  that 
apart  from  the  uncertainty  coming  from  the  assumed  capture  probability  this  model  has 
no  free  parameters  and  that  the  result  agrees  nicely  with  the  experiment. 

In  order  to  substantiate  this  we  performed  classical  3D  Molecular  Dynamics  calcu- 
lations with  the  methods  of  refs.  [8], [9]  in  the  geometry  of  an  infinite  straight  cylinder. 
Here  30  particles  are  distributed  in  a disk  over  two  strings  with  different  occupation. 
Their  initial  positions  and  velocities  are  distributed  randomly  according  to  their  avail- 
able longitudinal  and  transverse  kinetic  energies  (the  number  of  particles  in  the  simula- 
tion is  limited  by  the  computer  memory  for  the  fast  Fourier  transform).  The  small  disk 
of  the  cylinder  has  repeating  boundary  conditions,  thus  summing  up  all  Coulomb  inter- 
actions by  Ewald  summation,  see  ref.  [8].  The  effect  of  different  masses  or  frequencies 
in  this  straight  geometry  is  transformed  to  different  velocities  of  the  strings  according  to 
eq.  1,  i.e.  a relative  velocity  in  the  rest  frame  which  corresponds  to  a frequency  differ- 
ence of  75  Hz.  No  cooling  mechanism  is  incorporated.  All  positions  over  5000  betatron 
oscillations  in  a lattice  with  solenoidal  focusing  with  the  mean  tune  of  the  ESR  are  then 
Fourier  transformed  and  each  spectrum  has  been  averaged  over  1 1 runs  with  different 
random  initial  conditions  to  yield  Fig.  3.  Here  the  occupation  numbers  of  the  strings  are 
written  at  the  bases  and  the  relative  temperatures  with  respect  to  the  average  Coulomb 
interaction  energy  between  two  ions  are  indicated  in  the  key. 

The  upper  frame  treats  the  unlocked  case  where  the  distance  between  the  strings  is 
much  larger  than  the  thermal  diameter.  Here  the  difference  of  75  Hz  between  the  fre- 
quencies of  the  two  central  peaks  in  the  lower  curve  with  equal  occupation  corresponds 
to  the  relative  velocity  between  the  strings.  Although  the  strengths  of  the  peaks  vary  with 
the  occupation  their  distance  remains  constant.  In  the  bottom  frame,  the  locked  case,  the 
particles  have  lost  their  identity  resulting  in  a broad  random  distribution.  The  partially 
locked  case  is  intermediate. 

There  is  indication  that  the  effect  of  locking  or  capture  is  not  restricted  to  the  existence 
of  ultra-cold  chains.  As  pointed  out  in  ref.  [10],  during  the  preparation  of  the  precision 
mass  experiments  similar  anomalous  frequency  shifts  have  been  observed  in  an  experi- 
ment with  a bare  ’^^Au^^"*"  beam  at  295  MeV/u.  Its  isotopes  196,195, i94^i^79-(-  other 
strongly  populated  fragments  had  velocity  spreads  of  about  5v/v  = 10“^  and  thermal 
widths  of  about  4 mm.  By  scraping  the  primary  beam  down  to  a a width  of  600  fim  and 
under  steady  electron  cooling  the  velocity  spread  decreased  by  one  order  of  magnitude. 
Nearby  fragment  nuclides  with  slightly  larger  masses  then  acquired  anomalous  nega- 
tive frequency  shifts  of  the  order  of  a few  Hertz,  see  Fig.  4.  This  effect  remained  if  the 
primary  beam  was  further  scraped  to  a width  of  50  fim. 

Evidently,  the  low  intensity  fragment  isotopes  (with  negative  shifts  and  error  bars) 
and  masses  slightly  larger  than  the  mass  of  the  primary  beam  (with  nearly  zero  shifts 
and  no  error  bars)  have  been  absorbed  into  the  cloud  of  the  primary  beam  itself,  thereby 
reducing  their  frequencies  slightly.  Unfortunately,  there  are  no  isotopes  with  masses 
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FIGURE  4.  Measured  anomalous  frequency  shifts  in  a relatively  warm  beam.  The  dots  with  nearly 
zero  shifts  belong  to  isotopes  of  the  prmary  beam,  others  are  nearby  fragment  nuclides.  The  gray  dots  are 
questionable  by  the  lack  of  intensity  (after  ref.[  10]) 


slightly  less  than  the  primary  beam  masses  with  sufficient  population.  Hence  in  Fig.  4 
the  gray  dots  should  be  discarded  because  the  signals  came  from  very  few  ions  only. 

In  summary,  we  have  shown  that  with  the  high  precision  mass  measurements  in  the 
ESR  very  cold  multiple  strings  of  ions  have  been  produced  which  either  run  parallel  to 
each  other  or  lock  into  a single  string  if  their  displacement  becomes  very  small. 

The  Windows  simulation  program  RODS  (Representation  Of  Dynamical  Strings)  is 
available  on  the  web  page  wvw-aix. gsi  . de/~hasse/hasse . html.  The  author 
likes  to  thank  Z.  Patyk  for  pointing  out  this  problem  and  H.  Geissel,  K.E.G.  Lobner,  T. 
Radon  and  B.  Schlitt  for  valuable  discussions  on  the  experiments. 
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